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Abstract 



We present new Herschel/PACS images at 70, 100, and 160yum of the well-known, nearby, carbon-rich asymptotic giant branch star 
IRC+ 10216 revealing multiple dust shells in its circumstellar envelope. For the first time, dust shells (or arcs) are detected until 320''. 
The almost spherical shells are non-concentric and have an angular extent between ~40° and ~200°. The shells have a typical width of 
, and the shell separation varies in the range of ~10''-35'', corresponding to ~500-l 700 yr. Local density variations within one 
arc are visible. The shell/intershell density contrast is typically ~4, and the arcs contain some 50% more dust mass than the smooth 
envelope. The observed (nested) arcs record the mass-loss history over the past 16 000 yr, but Rayleigh-Taylor and Kelvin-Helmholtz 
instabilities in the turbulent astropause and astrosheath will erase any signature of the mass-loss history for at least the first 200 000 yr 
of mass loss. Accounting for the bowshock structure, the envelope mass around IRC+ 10216 contains >2Mo of gas and dust mass. It 
is argued that the origin of the shells is related to non-isotropic mass-loss events and clumpy dust formation. 

Key words. Stars: AGB and post-AGB, Stars: mass loss. Stars: circumstellar matter. Stars: carbon. Stars: individual: IRC +10216 



1. Introduction 

Stars of low and intermediate initial mass lose a significant frac- 
tion of their mass when evolving along the asymptotic giant 
branch (AGB). Pulsations in combination with radiation pres- 
sure on newly formed dust grains are thought to be chiefly 
responsible for driving the stellar wind. During the mass-loss 
event, a large circumstellar envelope (CSE) is formed. The mass- 
loss history of the stars is imprinted, in part, into the structure of 
the extended CSE. In the past decade, more and more evidence 
has been found that the envelope of many AGBs show deviations 
from spherical symmetry and/or exhibit temporal density var ia- 
tions (iNeri et al .l ll998l: ISchoier et al.ll2005l: iDecin et al.l[2QQ7h . 

IRC +1 0216 (CW Leo) is t he nearest (carbon-rich) AGB star 
(D~150pc,|C rosas & Menten|[T997b and serves as an archetype 
for the study of AGB stars. It has been the topic of many 
studies focussing on the structure of its thick envelope, which 
is created by a stellar \ yind carrying away the material at a 
rate of -1x10"^ Mo/yr (iDecin et al.l 1201 Oaf) . On smaU spatial 
scales (arcsec and below), the central regions of IRC +10216 
show the presence of a peanut- or bipolar- like structure (e.g., 
iMauron & Hugginsll2 QQQ: Menut et al. 2007^) together with a se- 
ries of clumps (e.g.,[Weigelt et al. 2002). Despite the presence 
of asy mmetric structure close to the star, iMauron & HugginsI 
) have detected multiple, almost concentric, shells (or arcs) 
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* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
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out to a distance of 55'' in radius on top of the smooth ex- 
tended envelope. The observed arcs are thought to be the re- 
sult of ambient Galactic light scattered by the dusty envelope, 
of which the shell/intersh ell density contrast is about of three 
(iMauron & Hugginsll200Ql) . The shells are incomplete and cover 
-30° -90° in azimuth with respect to the central star. Deep V- 
band images taken by iLeao et aP (l2006h show that the nested 
shells appear to extend to ~80'' and are composed of thinner 
elongated shells. The shells have a typical width of -1.6'' and 
are separated by ~5"-20", which corresponds to -250-1 000 yr. 

In this paper, we present new infrared images (at 70, 100, and 
160yum) of the dusty envelope around IRC+ 10216. Eor the first 
time, multiple dust shells are detected up to a distance of 320" 
from the central star. Although limited by the sensitivity of the 
data, these observations suggest that the shells are p resent all the 
way t o the bowshock at -450" at lOOyum (see also lLadjal et al.l 
I201QI) . In Sect. [2l we present the Herschel/PACS observations. 
We analyse, in Sect. [3l the morphology of the dusty envelope 
around IRC+10216 and the compatibility of the PACS infrared 
observations with optical data. We then discuss, in Sect. IH the 
density variations in the arcs, the mass-loss history and the pos- 
sible origin of the multiple shells. Einally, the conclusions are 
presented in Sect. \5\ 



2. Observations and data reduction 



The observations were obtained usin g the Photodetector A rray 
Camera and Spectrometer (PACS. IPoglitsch et aP |201Q|) on 
board the Herschel satellite (|Pilbrattetal.ll201Ql) . and are part of 



1 



L. Decin et aL: Multiple dust shells beyond 1 arcmin in the CSE of IRC+ 10216 



NOT DECONVOLVED DECONVOLVED 




Figure 1. PACS images of IRC + 10216 at 70 jim (upper), 100 fim (middle), and 160 fim (bottom) after subtraction of the smooth halo 
of the circumstellar envelope. Left-hand column shows the data without deconvolution, right-hand column when deconvolution is 
applied before the subtraction of the smooth envelope structure. The field-of-view (FOY) for the left-hand column is 16'x 11', while 
for the deconvolved images the FOY is lO'x 10'. North is at the top, east is to the left. The pixel size is V for the 70 and lOOyum 
images and is for the 160yum image. The six radial spikes are due to the support spider structure of the secondary mirror; these 
regions are not taken into account in our deeper analysis (see also Fig. [4]). The diagonal striping in the 100 jim in iage is instrumental . 
The innermost arc is found at 24'', the outermost one at 320". The bowshock detected at 480" to the east by lLadjal et al.l (l2010h 
can be seen in the left-hand image at 160yum. 
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the M ESS guaranteed time key programme (iGroenewegen et all 
l20TTh . which is investigating the dust and gas chemistry of a 
large sample of post-main-sequence stars. The scan-map observ- 
ing mode was used with a scan speed of 20'7s. The total sky 
coverage is 30'x30'. Two scan maps were taken with a scan- 
ning angle of 90° between the two scans to achieve a homoge- 
neous coverage. The P ACS 160yum data were already presented 
by iLadjal et aD (|201Q|) , who discussed the detection of a bow- 
shock at -585'' when measured from the nearest focus of an el- 
lipse fit to the shape of the bowshock, or at 480'' when measured 
from the central star. 

The data reduction was performed using the 
Herschel Interacti v e Pro cessing Environment (HIPE) (see 
IGroenewegen et al.] l201ll for more details). To remove 

low-frequency noise, causing additive brightness drifts, 

the scANAMORPHOs routinc (Roussel 1201 ll) was applied. 
Deconvoluti on of the data was d one using the method as 
described in lOttensamer et al.l (1201 ll) . 

To emphasize the shell morphology, we removed the central 
extended envelope halo by subtracting a smooth, azimuthally av- 
eraged profile represented by a power law r"^ for r < 500". The 
power a was determined directly from the local minima in the 
data for each image separately, and represents the azimuthally 
averaged (pseudo-) continuum intensity. The power index a 
varies between 2.07 and 2.50. The power index is diff'erent from 
that for a simplified r"^ power law, obtained in the case of an 
optically thin, homogeneous, spherical envelope with a constant 
dust temperature or a r"^^ power law when the dust temperature 
is proportional to r"^-^ . As discussed in more detail in Sect. 14.11 
the smooth underlying intensity represents both the continuum 
intensity (see red dotted line in Fig. [S]) and the additive contri- 
bution from each shell, resulting in a pseudo-continuum with a 
higher value for a. 

While the PACS instrument off'ers a pixel size of 3.2" in the 
70 and 100 jim bands and of 6.4" in the 160 yum band, the images 
shown in this paper are oversampled by a factor of 3.2, resulting 
in a sampling of 1" and 2" per pixel, respectively. The PACS 
resolution at 70, 100, and 160yum is, respectively, ~5.5", -6.7", 
and -11". The final PACS images are shown in Fig.[T] 

3. Results 

It is clear that the extended envelope of IRC-i-10216 shown in 
Fig. [T] is composed of multiple shells (or arcs). Thanks to the 
high sensitivity of PACS, at least a dozen shells are detected for 
the first time beyond a radius of 1 arcmin, with the outermost 
shell being located at -320" fr om the central star For a terminal 
gas vel ocitv. v.^.. of 14.5 km/s ([Deci n et al. 2010bV a distance of 
150pc (ICrosas & Mentenll 1997b . and neglecting the drift veloc- 
ity (see Sect. 14.3b , the outermost arc was ejected about 16 OOOyr 
ago. The shell separation is typically between -10" and 35" (or 
500-1 700 yr), but some shells intersect each other. Owing to the 
complex point- spread function (PSF) and the bright contribution 
from the central target at these infrared wavelengths, the inner- 
most shell detecte d by PACS is (only) at -24"(see Fig.[2l), while 
iLeao et aP (l2006l) detected the first shells as close as 4" from 
the central target using V-band images. A good correspondence 
is seen between the innermost arcs seen in the PA CS images 
and th e outermost shells seen in the V-band data of iLeao et al.l 
(l20Q6l) (see Fig. O. As was the case for the i nner shells (<1', 
e.g. lMauron & Hugginsl[T999l:lLeao et al.ll2006l) . all shells in the 
PACS images appear to be non-concentric and azimuthally in- 
complete. Some arcs have a positive and others a negative incli- 
nation. However, while the arcs seen in the V-band image cover 



-30°-90° in azimuth, the arcs in the PACS 100 yum image extend 
-45°-200°. The source IRC-hl0216 was also observed with 
the SPIRE instrument in the framework of the MESS GTKR 
However, after applying the same reduction strategy, only two 
faint arcs at -1' and 2' are visible in the SPIRE 250 yum image 
(see Fig. lA.ll in the Appendix). 




Figure 2. Zoom into the inner 5'30" x 4'30" of the PACS 
lOOyum image (cf. the middle panel in the right-hand column 
of Fig.[T]). The red circle with a radius of 15" marks the region 
where some PSF artefacts are still visible. 



The arcs seen i n the B and V band images p resented by 
iMauron & HugginsI (119991) and iLeao et al.l (l2006l) are caused 
by dust scattering of the external illumination by the ambient 
Galactic light, but the PACS images show the direct emission of 
dust grains of which the main heating source is the thermal emis- 
sion of nearby dust grains, with minor contributions from either 
the stellar UV flux (close to the star) or the diff'use Galactic light 
(in the outer envelope). The highest contrast between the arcs 
is reached in the PACS lOOyum image. Applying aperture pho- 
tometry on the azimuthally averaged deconvolved PACS images 
and using a modified blackbody of the form By • A~l^ , as expected 
for a grain emissivity Qabs ~ ^"^ with B equa l to 1.2 (represent- 
ing amorphous carbon, lMennella et al.ll 19981) . we derived a dust 
temperature between 108+5 K at 20" and 40+5 K at 180". The 
tempe rature of the bowshock at 480" is 25 + 3 K (iLadjal et al.l 
120101) . 

The thickness of the shells typically varies within the 
range -5"-8", and most of them are hence not resolved. 
Local variations in one arc are se en. The typ i cal thi ckness de- 
rived from the V-band image bv iLeao et al] (l2006l) is around 
1.6" +0.3 "(obtained for a pixel size of 0.2" and a spatial res- 
olution of -0.6"), and no relation was found between the shell 
thickness and the distance to the center, which could be expected 
from the natural expansion of the envelope. 

iDinh-V-Trung & LimI (l2008l) found that the gas shells seem 
to be more concentrated on the western part of the envelope and 
scarce at position angles of -0°-30° (measured counterclock- 
wise from north). The PACS observations confirm the higher 
concentration of shells in the western part, but the lowest num- 
ber of arcs close to the target are found in the eastern part (and 
not in the north-eastern region; see also Fig. [5]). 
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Figure 3. Comparison between PACS images and deconvolved VLT/FORSl V-band image (see Fig. 4 in iLeao et al.ll2006l) . The 
field-of-view is 204''x 204''. Left-hand i ma^e: Deconvolve d PACS lOOyum image (flux units being mJy/arcsec^) with superposed 
contours (in blue) of the V-band image by lLeao eTal ] (l2006h . The red circle marks the region in which the flux is unreliable because 
of PSF artefacts (see Fig. O. Right-hand image: Combined image of the PACS 70yum (green), PACS 100 yum (red), and V-band 
image (blue). The outermost arcs detected in the V-band image (at a radius of -32'', 53", and 73") are clearly seen in the PACS 
images as well. (For the best result, viewing this figure electronically is recommended). 



4. Discussion 

4.1. Shell intensity and density contrast of the WOfim image 

From the PACS images, the relative intensity variations in the 
arcs can be derived. This is illustrated in the case of the lOOyum 
image where the arcs are most clearly visible. We only used re- 
gions unaff'ected by PSF features caused by the support spider 
structure of the secondary mirror (see Fig. |4] and Table [T]). To 
measure the intensity variations along the radial intensity pro- 
files between the diff'erent local minima and maxima, the PACS 
images were p rojected from pola r to cartesian coordinates (see, 
e.g.. Fig. 6 in iLeao et al ] l2006h . The mean intensity enhance- 
ments in the difl'erent sectors are displayed in Fig. [3 The ratio 
of the integrated intensity to the full envelope intensity in the 
individual shells varies between ~2 and 15%, with the brightest 
arcs seen in sector 0. The lowest numbers of arcs is seen in sec- 
tor 5. While the typical thickness of the shells varies within the 
range ~5"-8", the shells appear broader in Fig.[5]because of the 
non-concentric nature of the shells. 



Table 1. Position angles (PA; measured counterclockwise from 
north) of the difl'erent sectors used in the analysis of the intensity 
variations in the arcs (see Fig. [4]). 



Name 


PA 


Name 


PA 


Sector 


320°-330° 


Sector 4 


125°-150° 


Sector 1 


300°-330° 


Sector 5 


65°-95° 


Sector 2 


255°-280° 


Sector 6 


17.5°-35° 


Sector 3 


185°-215° 


Sector 7 


0°-17.5° 




Figure 4. PACS lOOyum image (in units of mJy/arcsec^) of 
IRC -hi 02 16 after subtraction of the smooth halo of the circum- 
stellar envelope. The field-of-view is 975 "x 660". North is at the 
top, east is to the left. The pixel size is 1". The central target is 
marked by the innermost white circle. The six radial spikes indi- 
cated by the white sectors are due to the support spider structure 
of the secondary mirror; these regions are not taken into account 
in the analysis. The position angles of each sector used in the 
analysis are listed in Table [T] 



As shown by iMauron & HugginsI (l2000h . the shell density 
contrast can be obtained from the shell intensity profiles. To 
illustrate this, we modelled the extended radial intensity pro- 
file taken at a position angle (PA) of 325° in the 100 yum image 
(see Fig.O. The smooth envelope structure is modelled assum- 
ing a constant mass-loss rate, i.e. the density varies as r"^ (or 
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Figure 5. Mean intensity enhancements (in Jy/arcsec^) in each 
sector of the PACS 100 yum image (see Fig.|4|. The position angle 
covered by each sector is listed in Table [T] The different sectors 
are offset to each other. The upper panel shows the intensity (in 
Jy/arcsec^) in the arcs between 65'' and 350'' away from the cen- 
tral target, while the lower panel shows the outermost ar c s unti l 
600'' away. The bowshock as discussed by iLadjal et aP (l2010l) 
is seen at -450" in sector 5 (see also Sect.|B]in the Appendix). 



p(r) = a • r"^, with a = M/(47rVco)). The shell-intershell contrast 
is characterized by p(r) = f(r) - a - r~^, where /(r) represents 
the relative change i n den sity in each region. As discussed by 
iMauron & HugginsI (l2000h . the modelling assumption that the 
shells are complete, might overestimate their contribution to the 
radial profile since the arcs tend to be prominent over a limited 
azimuthal angle. This is, however, roughly compensated by the 
additional contribution of the other shells which are not seen as 
prominent arcs when they do not intersect the plane of the sky. 
The thermal emission of the dust grains is modelled assuming a 
spherically symmetric envelope and LTE (local thermodynamic 
equilibrium), i.e. the intensity ly can be written as 



+ Sy(l-e-^^), 



with /bg being t he interstellar radiation field taken from 
iPorter et aP (l2006l) . Ty the optical depth along the line of sight, 
and Sy the source function assumed to be characterized by 
^y(^dust(^))- The fractional contribution of the first term on the 
right-hand side of Eq. ([T]) to the second term is <0.4% for 
r < 350", and can therefore be neglected. The dust tempera- 



ture distribution is taken from the SED modelling presented in 
iDecin et aP (l2010ah . and can be approximated as 



Tdmt(r) = 848 ^ (r/rinner)' 



-0.48 



(2) 



with Tinner the dust Condensation radius at 0.127". This smooth 
(analytical) dust temperature distribution yields values lower 
than the dust temperatures derived in the shells from the aper- 
ture photometry in Sect. 3. This is because the functional form 
of a modified blackbody is only valid in the case of optically 
thin emission assuming a constant dust temperature along the 
li ne-of- sight. The mor e sophisticated dust modelling presented 
in iDecin et aP (1201 Oah takes optical depth effects into account 
when considering a radially varying dust temperature. 
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Figure 6. Multiple- shell radiative-transfer model fit to the radial 
intensity profile at a PA of 325° in the 100 yum image. Bottom 
panel: The black line shows the deconvolved observations, the 
red dotted line shows the predicted continuum intensity based 
on Eq. ([T]), and the green full line shows the multiple shell fit to 
the data based on Eq. ([T]) using the fractional density variations 
shown in the upper panel. Upper panel: Relative density change 
f(r) in each arc. The profile was not fitted for r < 15", because 
of the complex PSF structure in this region. 



The radial intensity profile shown in the bottom panel of 
Fig. [6l was fitted using the fractional density variation displayed 
in the upper panel of Fig.O Density contrasts of up to a factor of 
ten are seen in some arcs, but a more typical value is about four. 
This is in good agreement with the results deduc ed from the in- 
ner ar cs in the scattered V-band image by Maur on & HugginsI 
(l200Qh . Our results indicate that the arcs contain -50% 
dust mass than the smooth envelope. 



more 



4.2. Tracing the mass-loss history of IRC+10216 



Sahai & ChronopoulosI (l2010h . iLadial et al.l (l2010h . and 
Cox et al.l (12011 ) discussed the detection of a bowshock east- 
ward of the central star, which is the result of the interaction 
between the AGB wind and the inter stellar medium (ISM). 
From the ultraviolet (UV) images, ISahai & ChronopoulosI 
(I2OIOI) deduced the approximate locations of the termination 
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shocl^3, astropaus^l, and bowshocljl eastward of the central 
target, to be at -500'', 590'' and 670", respectively. 

The PACS images show a clear intensity enhancement 
around 450" at 100 yum (see Fig.© and aro und 480" at 160iUm 
(see Fig. lB.ll in the Appendix). As shown by lLadjal et al.l (l2010l) . 
a dust temperature of 25 ± 3 K is derived from the PACS and 
SPIRE images for a mixture of carbon-rich (from the stellar 
wind) and oxygen-rich (from the ISM) material (J3= 1.6). This 
suggests that the PACS images detect the turbulent astropause, 
where instabilities mix the ISM and circumstellar material and 
that the termination shock is situated < 450". Hydrodynamica l 
simulations with the AMRVAC-code (Ivan Marie et al.l l201ll) 
show that the termination shock is indeed closer to the star, i.e. 
at 390" and that the bowshock is located at -650" (in accor- 
dance with the observations). The region interior to the termi- 
nation shock consists of the unshocked, freely streaming stellar 
wind. In the turbulent astropause, Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities dominate the picture. 

Although the very faint emission makes it impossible to de- 
tect dust shells beyond 320" in the PACS images, it is reasonable 
to assume that the nested shell structure extends all the way out 
to the bowshock (see Fig.O. The imprint of this complex mass- 
loss history will be traceable until the termination shock (at 390" 
or during a period of 19000yr). This part of the envelope has a 
mass of ~O.19M0 for a mean mass-loss rate of 1 x 10"^ Mo/yr. 
However, in the turbulent astropause region, any density con- 
trast between the shell and intershell regions will be wiped out 
by the instabilities. For a mean velocity of 1 km/s in the region 
between termination shock and bowshock (as deduced from the 
hydrodynamical simulations), this means that at least the first 
-200 000 yr of mass-loss history are lost. After undergoing mass 
loss for at least 219 000 yr, the total CSE mass is estimated to be 
higher than 2 M©, which is in accordance with the initial mass of 
3 < M < 5 Mo as estimated by Guelin et al. (1995). 

Using aperture photometry, one can also estimate the total 
dust (and gas) mass in the envelope, and compare that to the 
value of 0. 19 M© derived in the previous paragraph. The total — 
background- subtracted — flux in an aperture of radius 390" is 
-1750 Jy, of which some 30% comes from the radius <15". A 
first order estimate of the dust mass is 



(O l2005h with Fy the observed flux (in erg s"^ cm"^ Hz"^), A 
the wavelength (in cm), d the distance (cm), k the Boltzman con- 
stant (= 1.381 10"^^ erg K"^), T^ust the dust temperature (in K) 
taken to be 30 K (see Sect. [3]), and k^. the dust opacity at the 
wavelength obs erved. We adopted a dust absorptivity at lOOyum 
of 25 cm^ g"^ (lHildebrandlll983l), but other higher va lues are 
also found in the literature (e.g.. iMennella et aD 1 19981) . Using 
Eq. (O, the total estimated dust mass for the envelope within 
390 " is Mdust ~ 9 X l O'^^McT). Adopting a gas-to-dust ratio of 
250 (iDecin et al.ll2010ah . the total (gas-hdust) envelope mass (for 
r <390") is -0.23 Mq. Although in good agreement with the 
value of 0.19Mo derived in previous paragraph, one should re- 
alize that the uncertainty in this value is at least a factor of two 
owing to (1.) the complex PSF yielding an uncertain flux esti- 



^ where the stellar wind velocity goes from supersonic to subsonic 
values 

^ where the ram pressure of the ISM equals the ram pressure of the 
stellar wind 

^ where the ISM velocity goes from supersonic to subsonic values 



mate for r <15", and (2.) the uncertainties in dust opacities that 
depend strongly on the shape and the size of the grains. 

4.3. Origin of the sliells 

We note that similar circular shells are seen around t he carbon- 
rich post-AGB object AFGL 2688 (ISahai et al.ll 1998b . As noted 
already by other authors on the topic of the shell structure in 
IRC-hl0216, it is difl&cult to identify a physical mechanism that 
can explain the origin of the non-concentric incomplete shells, 
which occur at time separations between a few hundred yr and 
-2000 yr. This time interval is muc h longer than the normal stel- 
lar pulsation period of 649 days (iLe Bertrdll992h . and at the 
same time much shorter than the interpulse period between he- 
lium flashes, which for IRC-hl0216 is estimated to be between 
6 000 and 33 000 yr (Ladjal et al. 2010). 

Three difl'erent types of models have been proposed in the lit- 
erature. One type is based on the e ffects of a binary companion . 
However, as discussed in detail by iMauron & HugginsI (l2000l) . 
the irregular arc spacing seen in the optical and infrared images 
does not support this type of model. 

A different scenario is proposed by ISoked (l2000l) . who sug- 
gested that a solar-like magnetic activity cycle is behind the for- 
mation of the multiple arcs. Owing to the appearance of mag- 
netic cool spots during the active phase, the gas temperature 
can be locally reduced. This can enhance dust formation and 
thus lead to a higher mass-loss rate from magnetic cool spots. 
The magnetic cool spots cover only a fraction of the stellar 
surface, and this inhomogeneity might explain the dumpiness 
seen on smaller spatial scales and in the arcs. Starspots have 
not yet been observed on AGB stars, although magnetic field 
stren gths around only a few l a te-type stars have been m easured 
(^e.g. JVlem mings et al.|[2005b . ISoker & ClavtonI (Il999l) suggest 
that cool magnetic spots on the surface of AGB stars might have 
a radius of <0.1R^. This would imply that several large starspots 
should co-exist at the same time s ufficiently close to each other 
to explain the result obtained by iDinh-V-Trung & LimI (l2008l) 
that each shell typically covers 10% of the stellar surface at the 
time of ejection. However, taking dynamical effects into account, 
this starspot model would lead to a periodicity of the shells of 
half of the magnetic cycle (iGarcia-Segura et aDl200ll) . This pe- 
rio dicity contradicts th e irr egular shell spacings obs erved. 

ISimis et al.l (l200lh and lSandin & Hofned (l2004b studied the 
effect of drift-dependent dust formation. An intricate nonlin- 
ear interplay between gas-grain drift, grain nucleation, radia- 
tion pressure, and envelope hydrodynamics can result in gas 
and dust density shells that oc cur at irregii l ar per iods of a few 
hundred years. As explained by lSimis et al ] (l2001h . the shells in 
their models appear perfectly round owing to the assumption of 
spherical symmetry. The fact that the chemical-dynamical sys- 
tem regulating the behaviour of the envelope is extremely stiff 
and reacts violently to all kind of changes suggests that in nature 
the shells will be incomplete or clumpy. However, the predicted 
shell density variations are much larger than the observed ones. 
This might be because dumpiness is not accounted for in this 
spherically symmetric models (see also next paragraph). 

Some constraints are set by the (observational) link between 
the inner regions, which for both IRC -hl0216 and AFGL 2688 
display a peanut- or bipolar-like structure in a clumpy environ- 
ment and the spherical-like shel ls observed on larger scales. 
According to iMenut etaP (l2007l) . the change from roughness 
on small scales to smoothness on larger scales can be explained 
by the dilution of the dust clumps while they move outward, 
or a shadowing effect in combination with multi- scattering of 
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the stellar and dust-emitted photons. In addition, backwarming 
effects during non-isotropic dust formation will inhibit the for- 
mation of new dust clumps between the star and already formed 
dust clumps. However, when a dust clump has moved outward or 
in any other direction where the number of nearby dust clumps is 
lower, dust formation is still possible. Repeated over a longer pe- 
riod, this process might result in spherical-like shells. This idea 
is supported by the high angular resolution observations of the 
HC3N and HC5N mo lecules, which have a cluni py and quasi- 
spherical distribution (iDinh-V-Trung & L"imll2008l) . The location 
of the arcs of the cyanopolyyne molecules coincide very clo sely 
with the du st arcs iden t ified b y iMauron & HugginsI (Il999l) and 
later on by iLeao et aP (l2006l) . suggesting that the drift veloc- 
ity is zero during thes e phases of higher dus t and g as density. 
From the observations, iDinh-V-Trung & LimI (120081) concluded 
that the arcs are true three-dimensional expanding features, and 
not projections on the plane of the sky of just a few l arger three- 
dimensional shells as suggested by lFong et al.l (l2003l) . Moreover, 
their velocity channel maps suggest that there might be slight 
variations in either the ejection velocity or ejection time between 
diff'erent parts of the shell, which might explain the non-circular 
nature of the arcs. 

The R\CS observations in combination with already pub- 
lished optical, infrared, and sub-millimetre images of the enve- 
lope around IRC -hi 02 16 suggest again that the dust formation 
and the onset of the mass loss is highly anisotropic, and that 
local density enhancements prevail until the termination shock. 
Local thermodynamical and sometimes magnetic fluctuations 
might have a large impact on the in- situ formation of new gas 
and dust species. The formation of one clump will have an ef- 
fect on its close environment, more specifically on the nucle- 
ation and growth efficiency of new dust clumps. Dilution and 
multi- scattering eff'ect will reduce the original contrast bright- 
ness of the clumps, but the continuous formation of new clumps 
combined with the eff'ects of non-equilibrium drift-dependent 
dust formation will result in shell-like structures as observed on 
larger scales. However, even within these shells, homogeneity 
will not prevail but relics from the original dust and gas clumps 
will remain. This complex envelope structure and local den- 
sity enhancements should be incorporated into chemical mod- 
els that try to explain the formation of complex molecules in 
the envelope around IRC-hl0216 (e.g ICordiner & Milladl2QQ9l: 
lAgundez et ani2010l:ICherchneffll201 Ih . 

5. Conclusions 

We have presented Herschel R\CS observations at 70, 100, and 
160yum that detect for the first time dust shells (or arcs) out 
to a radial distance of 320'' in the envelope around the well- 
known carbon-rich AGB star IRC-hl0216. The shells are non- 
concentric and azimuthally incomplete. The shell/inter shell den- 
sity contrast is typically a factor of about four. Within an arc, 
local density variations are clearly visible. The arcs represent 
true density variations in the circumstellar envelope. From the 
R\CS image, we have been able to infer the mass-loss history 
during the past 16 OOOyr. It is argued that the nested shell struc- 
ture should at least be visible out to the termination shock, at 
-390''. Instabilities in the turbulent astropause would wipe out 
any density contrast. 

Acknowledgements. We are grateful to 1. C. Leao for providing the V-band 
image of IRC+ 10216. This research has made use of NASA's Astrophysics 
Data System Bibliographic Service and the SIMBAD database, operated at 
CDS, Strasbourg, France. PACS has been developed by a consortium of in- 
stitutes led by MPE (Germany) and including UVIE (Austria); KUL, CSL, 



IMEC (Belgium); CEA, OAMP (France); MPIA (Germany); IFSI, OAP/AOT, 
OAA/CAISMI, LENS, SISSA (Italy); TAG (Spain). This development has been 
supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), 
GEA/GNES (France), DLR (Germany), AST (Italy), and GIGT/MGT (Spain). LD 
acknowledges financial support from the Fund for Scientific Research - Flanders 
(FWO). PR, BV, JB, and MG acknowledge support from the Belgian Federal 
Science Policy Office via the PRODEX Programme of ESA. FK acknowledges 
funding by the Austrian Science Fund FWF under project number P23586-N16, 
RO under project number I163-N16. 

References 

Agundez, M., Gernicharo, J., & Guelin, M. 2010, ApJ, 724, L133 
Gherchneff, I. 2011, A&A, 526, Lll 
Gordiner, M. A. & Millar, T. J. 2009, ApJ, 697, 68 
Gox, N., Decin, L., Kerschbaum, F., et al. 2011, A&A, submitted 
Grosas, M. & Menten, K. M. 1997, ApJ, 483, 913 
Decin, L., Agundez, M., Barlow, M. J., et al. 2010a, Nature, 467, 64 
Decin, L., Gernicharo, J., Barlow, M. J., et al. 2010b, A&A, 518, L143 
Decin, L., Hony, S., de Koter, A., et al. 2007, A&A, 475, 233 
Dinh-V-Trung & Lim, J. 2008, ApJ, 678, 303 
Fong, D., Meixner, M., & Shah, R. Y. 2003, ApJ, 582, L39 
Garcia-Segura, G., Lopez, J. A., & Franco, J. 2001, ApJ, 560, 928 
Groenewegen, M. A. T., Waelkens, G., Barlow, M. J., et al. 2011, A&A, 526, 
A162 

Guelin, M., Forestini, M., Valiron, R, et al. 1995, A&A, 297, 183 
Hildebrand, R. H. 1983, QJRAS, 24, 267 

Ladjal, D., Barlow, M. J., Groenewegen, M. A. T., et al. 2010, A&A, 518, L141 
Leao, I. G., de Laverny, P., Mekarnia, D., de Medeiros, J. R., & Vandame, B. 

2006, A&A, 455, 187 
Le Bertre, T. 1992, A&AS, 94, 377 

Li, A. 2005, in American Institute of Physics Gonference Series, Vol. 761, The 

Spectral Energy Distributions of Gas-Rich Galaxies: Gonfronting Models 

with Data, ed. G. G. Popescu & R. J. Tuffs, 123-133 
Mauron, N. & Huggins, R J. 1999, A&A, 349, 203 
Mauron, N. & Huggins, R J. 2000, A&A, 359, 707 
Mennella, V., Brucato, J. R., GolangeH, L., et al. 1998, ApJ, 496, 1058 
Menut, J., Gendron, E., Schartmann, M., et al. 2007, MNRAS, 376, L6 
Neri, R., Kahane, G., Lucas, R., Bujarrabal, V., & Loup, G. 1998, A&AS, 130, 1 
Ottensamer, R., Luntzer, A., Mecina, M., et al. 2011, in Why Galaxies 

Gare About AGB Stars: Shining Examples and Gommon Inhabidants, ed. 

F. Kerschbaum, Th. Lebzelter, & R. F. Wing, Astronomical Society of the 

Pacific Gonference Series, in press 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, LI 
Poglitsch, A., Waelkens, G., Geis, N., et al. 2010, A&A, 518, L2 
Porter, T. A., Moskalenko, I. V., & Strong, A. W. 2006, ApJ, 648, L29 
Roussel, H. 2011, A&A, submitted 
Sahai, R. & Ghronopoulos, G. K. 2010, ApJ, 711, L53 
Sahai, R., Trauger, J. T., Watson, A. M., et al. 1998, ApJ, 493, 301 
Sandin, G. & Hofner, S. 2004, A&A, 413, 789 
Schoier, E L., Lindqvist, M., & Olofsson, H. 2005, A&A, 436, 633 
Simis, Y. J. W., Icke, V., & Dominik, G. 2001, A&A, 371, 205 
Soker, N. 2000, ApJ, 540, 436 
Soker, N. & Glayton, G. G. 1999, MNRAS, 307, 993 

van Marie, A. J., Meliani, Z., Keppens, R., & Decin, L. 2011, ApJ, 734, L26 
Vlemmings, W. H. T., van Langevelde, H. J., & Diamond, R J. 2005, A&A, 434, 
1029 

Weigelt, G., Balega, Y. Y, Blocker, T., et al. 2002, A&A, 392, 131 



Appendix A: SPIRE 250 yum image of IRC +1 021 6 

The SPIRE 250 yum was reduced using the same reduction tech- 
nique as described in Sect.[2l Only two faint arcs, at -T and ~2' 
are seen. 



Appendix B: The turbulent interaction region with 
the ISM 

The detection of the bowshock/astropause (see Sect. 14.21) 
around IRC+10216 using the Herschel R\CS 160yum and 
SPIRE 250 | u m, 35 yum, and 500 yum images was published by 
iLadial et all (l2010h . To this detection, we now add the PACS 
lOOyum image, while the bowshock is not visible in the 70yum 
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Figure A.l. SPIRE 250yum image of IRC +10216 after subtrac- 
tion of the smooth halo of the circumstellar envelope. The field- 
of-view is 975''x 66C. North is at the top, east is to the left. 



image (see Fig. [T]). The bowshock is also visible in Fig. \T\ in 
the lOOyum image situated closer to the central target than in 
the 160 yum image. This is visualised in Fig. lB.ll At 100 yum, the 
bowshock structure is at -450'' from the central targe t, but it 
is at ~ 480'' in the 160 yum image. From Fig. 2 in Ladj al et aP 
(I2OIOI) . it seems that this trend is also seen for the longer wave- 
lengths covered by SPIRE. This might indicate that smaller (and 
warmer) dust grains are responsible for the lOOyum emission 
than at the longer wavelengths, in agreement with the theoret- 
ical predictions of .van Marie et al.. (.201 li) . 




Figure B.l. Mean intensity enhancements in the PACS 100 

(upper curve) and 160yum (lower curve) images in sector 5 
(see Fig. [4]) , zoom ing into the bowshock region discussed by 
iLadjal et aP (l2010h . When measured from the central target, the 
bowshock is clearly detected at 480'' in the 160yum image, and 
the corresponding feature is seen at 450" in the lOOyum image. 
The feature around 370" in the 100 yum image is not an 'arc' , but 
more of a blob-like enhancement as can be seen in Fig.|4l 
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